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ABSTRACT 

Precise Doppler searches for extrasolar planets find a surfeit of planets with 
orbital periods of 3-4 days, and no planets with orbital periods less than 3 days. 
The circumstellar distance, i?o, where small grains in a protoplanetary disk reach 
sublimation temperature (~ 1500 K) corresponds to a period of ~ 6 days. Interior 
to R , turbulent accretion due to magneto-rotational instability may evacuate 
the disk center. We suggest that planets with orbital periods of 3-4 days are so 
common because migrating planets halt once this evacuated region contains the 
sites of their exterior 2:1 Lindblad resonances. 

Subject headings: astrobiology - - circumstellar matter - - planetary systems: 
formation — planetary systems: protoplanetary disks — stars: formation 

1. Introduction 

Analytic calculations (Goldreich & Tremaine 1980; Ward 1997a) and numerical simu- 
lations (Nelson et al. 2000; Kley et al. 2001) suggest that protoplanets in a protoplanetary 
disk migrate rapidly into the star they orbit — so rapidly that it is a wonder any planets sur- 
vive at all. Small protoplanets torque the disk at Lindblad and corotation resonances, and 
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the resulting back-torque can propel a planet into the star in a matter of 10 5 (M P /M^y 1 
years (Ward 1997b). Ward (1997b) has dubbed this conundrum the Shiva problem, after the 
Hindu god of destruction. Large protoplanets may open a gap in the disk via their resonant 
torques, and so become locked to the disk's viscous spreading, a process which may dump 
the most of the disk and planet onto the star within 10 7 years or less, depending on the disk 
viscosity. 

Figure 1 shows the distribution of the orbital periods of the innermost Doppler planet 
candidates, summarizing data from Marcy et al. (2002) 2 . These candidate planets all orbit 
stars with masses in the range 0.7-1.4 M & . Figure 1 suggests that whatever mechanism 
halted the migration of these planets operates best at an orbital period of ~ 3 days and 
ceases to operate at shorter periods. Of the 20 planets with periods less than 20 days, 8 
have periods in the range 3-4 days. No planet has a period less than 2.98 days. This trend 
appears to be real and not an artifact of observational selection; the primary precise-doppler 
surveys are complete for Jupiter mass planets out to a period of > 0.5 years (Butler et al. 
2001). 

Occasionally the interactions among two planets and a star can leave a planet trapped 
by stellar tides into a circular orbit at ~ 0.04 AU (Rasio & Ford 1996). Sometimes an 
accreting planet may overflow its Roch lobe, losing mass to the star, and this process may 
halt the planet's migration (Trilling et al. 1998). But these phenomena appear to be rare. 

Scenarios in which planets migrate by interacting with a disk of planetesimals without 
gas (Murray et al. 1998) provide a disk truncation radius where planets may gather: the 
radius where planetesimals become hot enough to sublimate. The problem with gas-less 
migration schemes is that substantially changing the orbit of a Jupiter-mass planet requires 
roughly a Jupiter mass of planetesimals. 

In contrast, optically-thick disks with more than a Jupiter mass of gas appear to be 
ubiquitous around young stars. Lin et al. (1996) have suggested that planet migration ends 
where the gaseous protoplanetary disk meets the stellar magnetosphere. We suggest an alter- 
native explanation for the pile-up of planets near the 3-day period: a gas disk truncated at a 
temperature of 1500 K by the onset of Magneto- Rotational instability (MRI; Chandrasekhar 
1961; Balbus & Hawley 1991). 



2 This website about extrasolar planets is available at http://cxoplanets.org/ 
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Fig. 1. — Histogram of the orbital periods of the extrasolar planet candidates detected by 
the Precision-Doppler technique. 
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2. The Theory 

A planet with orbital period p and semimajor axis a interacts with a disk at Lindblad 
resonance sites, located where the natural epicycle period of the gas in the disk (roughly the 
Keplerian orbital period) is p(m + l)/m. The most distant Lindblad resonances from the 
planet are the m — 1 and m = —2, located at roughly 1.59a and 0.63a respectively. The 
planet torques the disk at these resonances, and the back-torque from the inner Lindblad 
resonances tends to add angular momentum to the planet's orbit, while the back-torque 
from the outer Lindblad resonances tends to remove angular momentum to the planet's 
orbit (Goldreich & Tremaine 1980). However, over a wide range of disk pressure and density 
gradients, the torques from the outer Lindblad resonances dominate the torques from the 
inner Lindblad resonances, causing the planet to spiral into the star (Ward 1997a). 

Consider a planet contained in a protoplanetary disk with an inner truncation radius 
of Rq, corresponding to an orbital period of Pq. When the planet migrates inward to where 
it has an orbital period p < 2P ( a ^ 1.59i?o), it may begin to migrate inwards faster 
than normal because its interior Lindblad resonance sites lie in the evacuated regions, so 
the outward torques from these resonances disappear. When the planet reaches an orbital 
period of 1/2P (a = 0.63_R ), and all of the planet's inner and outer Lindblad resonances 
are contained in a central evacuated region, the disk-planet interaction may stall, ending 
the planet's migration (Lin et al. 1996). A planet with a gap is not immune; it must also 
stop migrating at a period of 1/2P , where it can no longer dynamically communicate with 
enough mass in disk material to substantially affect its orbit. 

We propose that protoplanetary disks are evacuated interior to radius i? because of the 
powerful magneto-rotational instability which afflicts conducting Keplerian disks. Proto- 
planetary disks may suddenly become conducting where they reach a temperature of 1500K 
and potassium ionizes (Stone et al. 2000) and dust sublimates, removing recomination sites 
(Sano, T. et al. 2000). Numerical simulations of saturated MRI show that it can provide a 
disk a viscosity (Lynden-Bell & Pringle 1974) in the range of a ~ 0.004 to a ~ 0.1 (Hawley 
et al. 1995; Ziegler & Riidiger 2001) depending on the assumed resistivity. Exterior to R , 
the MRI may be damped completely (Gammie 1996; Reyes-Ruiz 2001). 

In a steady state disk, the accretion rate M is uniform and continuity dictates that 
RYiVr = constant, where R is the radial coordinate, £ is the surface density, and vr is the 
radial drift velocity. In a viscous accretion disk, vr ~ v ■/ 'R, where v is the viscosity. So 
roughly speaking, £ oc i/ -1 , or using the a prescription for disk viscosity, £ oc a -1 . 

Magnetorotational instability can plausibly cause the effective a to climb by a factor of 
10-1000 interior to Rq, which translates into a comparable drop in surface density. When 
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planets reach a = 0.63P in the centers of disks evacuated by MRI-driven accretion, their 
migration will slow by this factor of 10-1000. The high viscosity in this zone also acts to 
inhibit gap formation and close up existing gaps around massive planets. 

In radiative-transfer models of protoplanetary disks, a surface layer of small grains 
intercepts the starlight first and rations it to the disk midplane (Calvet et al. 1991; Malbet 
& Bertout 1991; Chiang & Goldreich 1997; Sasselov & Lecar 2000). These small grains have 
hotter equilibrium temperatures than a blackbody at a given distance from the star, so they 
sublimate farther from the star than larger grains. We can find Rq by iteratively solving 
the radiative equilibrium equation for the temperature of the small grains as a function 
of distance from the star (Backman & Paresce 1993). We performed such a computation 
assuming ~ 0.1/xm grains — bodies whose emissive and absorptive efficiencies stay constant at 
wavelengths shorter than 0.1 /xm, but decline as A -1 at longer wavelengths. We find that the 
grains reach 1500 K at R = 0.067 AU (P = 6.3 days) for a solar type star, or R = 0.055 
AU (Po = 4.8 days) for star with solar luminosity and mass, but effective temperature of 
4000 K, a common model for a T Tauri star. 

Planet migration in a disk with such an inner truncation radius would halt at p ~ 
P /2 = 2.4 days; near the observed pile-up of precise-Doppler planets. The luminosity of 
a solar-mass star varies over roughly a factor of ~ 5 during the time from 1 million to 10 
million years after its birth (D'Antona & Mazzitelli 1994), so the value of Po relevant to 
halting the migration of a given planet could conceivably range over a factor of ~ 2 due to 
this effect. For comparison, if we assume a 1 Jupiter mass planet has a radius of 2.3 Jupiter 
radii at an age of 1 million years (Guillot & Showman 2002), we calculate that the orbital 
period where it overflows its Roche radius is ~ 1.3 days. 

3. Discussion 

3.1. OBSERVATIONS OF PROTOPLANETARY DISKS 

The likely hosts of protoplanetary disks are T Tauri stars (< 2M Q ) and Herbig Ae/Be 
stars (~ 2-8M Q ). Near-infrared photometry of Herbig Ae/Be stars reveals a curious pattern; 
many of these stars have two-humped infrared excesses-broad emission peaks in the mid- 
or far-infrared, and then second emission peaks near 3/im (Lada & Adams 1992). A model 
in which accretion luminosity powers the near-IR hump requires optically-thick gas interior 
to the dust-sublimation point, a feature not corroborated by the observed SEDs (Hartmann 
et al. 1993) or measurements of optical veiling (Ghandour et al. 1994). The near-IR humps 
appear to represent the thermal emission peak of a population of dust at a temperature of 



- 6- 



1500 K (Sylvester et al. 1996; Malfait et al. 1998). The prominence of these near-IR humps 
may indicate the importance of R in the structure of Herbig Ae disks. 

Further evidence for the structural importance of R comes from long-baseline optical 
interferometry. Millan-Gabet et al. (2001) resolved 11 Herbig Ae stars at baselines of 21 and 
38 m with the IOTA interferometer in the near- infrared. The visibilities of the 11 resolved 
targets and their near-IR excesses were well modeled by dusty rings with inner radii of 0.3-5 
AU, possibly corresponding to Rq. Measurements of the near-IR sizes of the central regions 
of other young stars, including some T Tauri stars, with the Palomar Testbed Interferometer 
(Akeson et al. 2000) and using a single Keck telescope with an aperture mask (Tuthill et al. 
2001; Danchi et al. 2001) support this association; the observed near-IR sizes show a trend 
of increasing with the luminosity of the young star, which is consistent with the increase of 
Ro with luminosity (J. Monnier & R. Millan Gabet, in preparation). 

Natta et al. (2001) and (Dullemond et al. 2001) have interpreted the photometric and 
interferometric observations of Herbig Ae/Be stars as signs that protoplanetary disks may 
have tall, passively-heated inner walls; hot dust on the wall provides the 3/xm emission peak. 
This inner wall could conceivably represent the zone where the opacity and the surface density 
suddenly drop in our model. Further high-resolution interferometric studies of Herbig Ae/Be 
and T Tauri disks should help clarify how disk structure affects planet migration. 

3.2. MAGNETIC DISK TRUNCATION AND THE X-WIND THEORY 

Shu et al. (2000) review a theory of the interiors of protoplanetary disks which centers 
on the interaction between the stellar magnetic field and the inner edge of the conducting 
gas disk. In this "X-wind" theory, the disk is replaced by magnetic accretion columns and a 
bipolar wind interior to the radius R x . Earlier explanations of halting planet migration (Lin 
et al. 1996) have pointed to R x as the critical radius rather than Rq. 

Directly measuring R x is difficult. However, in equilibrium, the stellar rotation period 
matches the Keplerian period at R x , so measuring the stellar rotation rate and assuming 
that T Tauri stars are rotationally locked to their disks may provide an estimate of R x . 
However, the process of breaking the stellar rotation is probably slow (Hartmann 2002), 
and the rotation periods of T Tauri stars range from 0.5 to > 8 days (Stassun et al. 1999). 
Consequently, it is hard to associate the 3-day period of the closest-in extrasolar planets 
with a special stellar rotation rate, or with R x . 

In the X-wind model of chondrule formation, a geometrically-thin, optically-thick disk 
of planetesimals interior to R x suffers mutual collisions and heating from magnetic flares 
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(Shu et al. 1996, 2001). Material from these planetesimals become proto-chondrules. In our 
picture, Type I migration helps deliver planetesimals to the magnetically-active zone. We 
imagine a ring of planetesimals collecting near 0.63i?o- Objects smaller than a few km may 
migrate in further than 0.63i? via plasma drag (Shu et al. 1997, 2001). 

3.3. TORQUE REVERSAL 

We briefly considered the following alternative mechanism for halting planet migration. 
Continuing the calculations of Ward (1997a) shows that with only a slight positive temper- 
ature gradient ((R/T)dT/dR > 0.6, i.e. the disk temperature increases with radius faster 
than roughly R° ), the torque on the planet from disk tides can reverse. Such a positive 
temperature gradient could conceivably occur where the disk regains sight of the star beyond 
the shadow of the disk's tall inner wall. However, the models of Dullemond et al. (2001) 
suggest that the positive temperature gradient in this region is largely confined to the surface 
layer of the disk, which does not have enough surface density itself to affect planet migration. 

4. Conclusion 

We have painted a picture of the central regions of protoplanetary disks in which a 
low-surface-density, rapidly accreting central region catches inwardly-migrating planets like 
flypaper. A good way to test our conjecture that disk temperature — not the stellar mag- 
netic field — determines the orbital radii of the innermost surviving planets is to extend the 
sample of planet search target stars to include a wider range of stellar masses. For example, 
we predict that planets around early A main sequence stars will collect at a radius much 
farther from the star (~ 0.3 AU) than planets around solar type stars. The high rotation 
rates of main sequence A stars betray that magnetic disk truncation can not rescue planets 
migrating into these stars until they reach near 0.03 AU. Doppler planet search techniques 
may not work for A stars because of their rotationally-broadened spectral lines. However, 
astrometric searches, transit searches, and direct-imaging searches for extrasolar planets are 
not so limited. 

We thank Joe Hahn, Doug Lin, Dimitar Sasselov, Hsien Shang and Bill Ward for helpful 
conversations. 
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